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In this paper we present the first analytic model for vorton forma- 
qq ■ tion. We start by deriving the microscopic string equations of motion in 

' Witten's superconducting model, and show that in the relevant chiral 

Q\ , limit these coincide with the ones obtained from the supersonic elastic 

models of Carter and Peter. We then numerically study a number of 
» I | solutions of these equations of motion and thereby suggest criteria for 

f> |- deciding whether a given superconducting loop configuration can form 

a vorton. Finally, using a recently developed model for the evolution of 
currents in superconducting strings we conjecture, by comparison with 
these criteria, that string networks formed at the GUT phase transition 
should produce no vortons. On the other hand, a network formed at 
the electroweak scale can produce vortons accounting for up to 6% of 
^ ■ the critical density. Some consequences of our results are discussed. 
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r"| ■ As first pointed out by Witten |5J], cosmic strings can in some circumstances 

Oh! (typically when the electromagnetic gauge invariance is broken inside the string) 

behave as 'superconducting wires' carrying large currents and charges — up to the 
| order of the string mass scale in appropriate units. The charge carriers can be 

. either bosons or fermions (see Q for a review). The former type occurs when it 

becomes energetically favourable for a charged Higgs field to have a non-zero vacuum 
expectation value in the string core; the latter happens when fermions couple to the 
■ string fields creating fermion zero modes. 

It is well known that arbitrarily large currents are not allowed — there is a critical 
value beyond which the current saturates. In other words, for large enough wind- 
ing number per unit length, the superconducting condensate is quenched down, 
suppressing the current flow. Also, the current can decay by magnetic flux-line tun- 
nelling; this can be used to impose constraints on allowed particle physics models. 

If superconducting strings carry currents, they must also carry charges of similar 
magnitude. This includes not only charges trapped at formation by the Kibble 
mechanism but also the ones due to string inter-commuting between regions of the 
string network with different currents. Just like with currents, charge densities 
cannot have arbitrarily large magnitude — there is a limit beyond which there will 
no longer be an energy barrier preventing the charge carriers from leaving the string. 

A rather important point is that the presence of charges on the string tends to 
counteract the current quenching effect discussed above. In fact, numerical simula- 
tions of contracting string loops at fixed charge and winding number have shown B 
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that a 'chiral' state with equal charge and current densities is approached as the loop 
contracts. In this limiting chiral case, quenching is in fact eliminated completely. 
This has several important consequences. Strings that have trapped charges as a 
consequence of a phase transition can become superconducting even if the formation 
of a condensate was otherwise energetically unfavoured. More importantly, a string 
with both a charge and a current density will have a non-zero angular momentum. 

In the cosmological context, these strings would of course interact with the cosmic 
plasma, originating a number of interesting consequences. The most remarkable of 
these, however, has to do with the evolution of string loops. If a superconducting 
string loop has an angular momentum, it is semi-classically conserved, and it tries 
to resist the loop's tension. This will at least increase the loop's lifetime. If the 
current is too large, charge carriers will leave the string accompanied by a burst 
of electromagnetic radiation, but otherwise it is possible that dynamically stable 
loops form. These are called vortons Q — they are stationary rings that do not ra- 
diate classically, and at large distances they look like point particles with quantised 
charge and angular momentum. Their cosmological significance comes from the 
fact that they provide very strong constraints on allowed particle physics models, 
since they behave like non-relativistic particles. According to current belief 
if they are formed at high enough energy scales they are as dangerous as magnetic 
monopoles, producing an over-density of matter in disagreement with observations. 
On the other hand, low-mass vortons could be a very interesting dark matter can- 
didate. Understanding the mechanisms behind formation and evolution is therefore 
an essential cosmological task. 

The overwhelming majority of the work done on cosmic strings so far was con- 
cerned with the structureless Goto-Nambu strings (but see and references therein 
for some exceptions). In the case of work on vortons, this means that somewhat 
ad-hoc estimates had to be made for some properties of the cosmic string network — 
notably for microscopic quantities such as current and charge densities. This is de- 
spite the fact it has been recognised a long time ago that, even though they might 
be computationally very useful Goto-Nambu models cannot realistically be 

expected to account for a number of cosmologically relevant phenomena, due to the 
very limited number of degrees of freedom available. Two such phenomena are the 
build-up of small-scale structure and charge and current densities. 

In this paper we fill this important gap by discussing the problem of vorton 
formation in the context of the superconducting string models of Witten and 
of Carter and Peter JlO| (sections [il] and III). Strangely enough, the issue of the 
conditions for vorton formation has been so far neglected with respect to those of 
their stability and cosmological consequences. We will start by introducing these 
models and determining the microscopic string equations of motion in each case. It 
will be shown that in the relevant chiral limit these equations coincide — this also 
provides the first conclusive evidence of the validity of the supersonic elastic models 
of Carter and Peter |I(J . 

We then proceed to study the evolution of a number of loop solutions of these 
equations numerically (sections IV and 0), and from the results of this analysis 
parameters will be introduced which characterise the loop's ability to evolve into a 
vorton state (section VI ) . Finally, we discuss a very simple phenomenological model 
for the evolution of the superconducting currents on the long cosmic string net work 



based on the dynamics of a 'superconducting correlation length' (sections VII 



|VIII ). Using this model we can therefore estimate the currents carried by string 
loops formed at all relevant times, and thus (in principle) decide if these can become 
vortons (section IX ) and calculate the corresponding density (section |x|) . 

Based on our results, we don't expect any GUT vortons to form at all. This is 
essentially because the friction-dominated epoch is very short for GUT-scale strings 
0, so their currents and charges are never large enough to prevent them from 
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becoming relativistic — and therefore liable to losses. Even if they did form, they 
wouldn't be in conflict with the standard cosmological scenario if they decayed soon 
after the end of the friction-domination epoch. 

Hence we conclude that, in contrast with previously existing estimates 
one cannot at the moment rule out GUT superconducting string models. We 
should point out at the outset that there are essentially three improvements in the 
present work which justify the different end result for GUT-scale strings. Firstly, by 
analysing simple (but physically relevant) loop solutions of the microscopic string 
equations of motion for the Witten model, we can get a much improved idea of how 
superconducting loops evolve and of how (and under which conditions) they reach 
a vorton state. Secondly, by using a simple model for the evolution of the currents 
on the long strings [ [Tl| we can accurately determine the typical currents on each 
string loop at the epoch of its formation. Finally, the use of the analytic formalism 
previously introduced by the present authors ^,|) allows us to use a quantitative 
description throughout the paper, and in particular to determine the loop sizes at 
formation. 

As will become clear below, when taken together these allow a detailed analysis 
of the process of vorton formation to be carried out, either in the Witten model (as 
is done in this paper) or any other that one considers relevant. In contrast, note 
that Davis & Shellard Q] restrict themselves to the particular case of the initial 
Brownian Vachaspati-Vilenkin loops with Kibble currents, and do not consider the 
subsequent evolution of the network. On the other hand, Brandenberger et al. Q 
make rather optimistic order-of-magnitude estimates about the process of relaxation 
into a vorton state. As it turns out, for high energy GUT scales, all these loops be- 
come relativistic before reaching a vorton state. Finally, neither of these treatments 
has the benefit of a quantitative model for the evolution of the long-string network 
|?|] which allows one to accurately describe the process of loop production. 

On the other hand, as we lower the string-forming energy scale we expect more 
and more efficient vorton production, and the 'old' scenario still holds. Therefore 
intermediate-scale superconducting strings are still ruled out, since they would lead 
to a universe becoming matter-dominated earlier than observationally allowed. Fi- 
nally, at low enough energy scales, vortons will be a dark matter candidate. For 
example, for a string network formed around T ~ 10 2 GeV (typical of the elec- 
troweak phase transition) they can provide up to 6% of the critical density. A more 
detailed discussion of these issues is left to a forthcoming publication [Q . 

Throughout this paper we will use fundamental units in which h = c = ks = 
Gm 2 Pl = 1. 



II. WITTEN'S MICROSCOPIC MODEL 

As first pointed out by Witten JlJ], a low-energy effective action for a supercon- 
ducting string can be derived in a way that is fairly similar to what is done in 
the Goto-Nambu case (see for example Q). One has to adopt the additional as- 
sumptions that the current is much smaller than the critical current and that the 
electromagnetic vector potential is slowly varying on the scale of the condensate 
thickness. 

The derivation then proceeds as in the neutral case, except for the use of the 
well-known fact that in two dimensions a conserved current can be written as the 
derivative of a scalar field. One obtains 



S = 



-Mo + \l° 




d 2 a 



(2.1) 
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/ d^xV^gF^F^ ; (2.2) 

the four terms are respectively the usual Goto-Nambu term, the inertia of the charge 
carriers, the current coupling to the electromagnetic potential and the external 
electromagnetic field (e ab is the alternating tensor); note that this applies to both 
the bosonic and the fermionic case 
Recalling the usual definitions 



Fab — FfivX^X ^ — A o a A a f) , 

and defining T a b to be the stress-energy tensor of the scalar field 
and the conserved current J a as 



(2.4) 



(2.5) 



J a = q- 



(2.6) 



we can obtain the following equations of motion by varying the action (2.2) with 
respect to A^, <fi and x M respectively 



F.^ = -Airqf = -4tt J d 2 ae ab x^ a (t>^S 4 (x - x(a a ) , 



(2.7) 



d a (V^TJ^b) + ^q^ ab F ab = , 



and 



d a 



-7 l a ° + -T a0 x\ b 



- 7 i ab + -r ab TZ p x° a x 
V Mo ' 

~(F a +n, 



where is the usual metric connection, 



(2.8) 

(2.9) 
(2.10) 



(2.H) 



and F a is the Lorentz force 



p a — F a T u T a 



We have also included the friction force 

f a = T 



„u> u, a J , 



(2.12) 



(2.13) 



using the same procedure as described in [[I3|j7j]. As shown in |ll|, plasma effects 
are subdominant, except possibly in the presence of background magnetic fields — 
either of 'primordial' origin or generated (typically by a dynamo mechanism) once 
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proto-galaxies have formed. Hence one expects Aharonov-Bohm scattering JT^J] to 
be the dominant effect, and consequently we have 

^ = wr (2 - 14) 

where Tj, is the background temperature and (3 is a numerical factor related to the 
number of particle species interacting with the string. 

The effect of self-inductance leads to the renormalisation of both the electromag- 
netic coupling and the scalar field cf>. Now, it is well known that the Maxwell-Faraday 
tensor includes both the external field and the field produced by the string itself, but 
it can be shown that if one follows this renormalisation procedure one can identify 
it with its external component, which henceforth we assume to vanish ||. 

As we already pointed out, when dealing with superconducting string loops we 
are essentially interested in the chiral limit of this model, that is 

7ab J a J b = 0, (2.15) 

or equivalently 

0' 2 = e 2 4> 2 . (2.16) 
In this limit, introducing the simplifying function <I> defined as 

*(X) = -^— , (2-17) 
Mo7oo 

and choosing the standard gauge conditions 

cr°=r, x-x' = 0, (2.18) 

(with dots and primes respectively denoting derivatives with respect to the time- 
like and space-like coordinates on the worldsheet as usual) the string equations of 
motion in an FRW background with the line element 

ds 2 = a 2 (dr 2 - dx 2 ) (2.19) 

(which implies that 700 = a 2 (l — x 2 )) have the form 

[e (1 + $)]'+ ^x 2 =$'-2-e$, (2.20) 



a 



and 



e(l + $)x+— (l-x 2 )x = 



(1-$) 



x'"" 



$ + 2-$ ) x' + 2$x' , (2.21) 
a ) 



where for simplicity we have introduced the 'damping length' 

l=a( 2 tf + I). (2.22) 

Finally, the worldsheet charge and current densities are respectively given by 

p w = qecf) , (2.23) 

and 

j w = q^- . (2.24) 

e 

Note that the Witten action is 'microscopic' in the sense of being built using 
only the properties of the underlying particle physics model [0 . In the next section 
we will analyse the equations of motion obtained form the action for the elastic 
supersonic models of Carter and Peter [HJ, which is is this sense 'macroscopic'. 
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III. SUPERSONIC ELASTIC MODELS 



In order to account for phenomena such as the build-up of charge and current 
densities on cosmic strings, one must introduce additional degrees of freedom on the 
string worldsheet. One such class of models, originally introduced by Carter and 
co-workers is usually referred to as elastic models (see M and references therein, on 
which the following two subsections are based). 



A. Basics of elastic models 



In general, elastic string models can be described by a Lagrangian density de- 
pending on the spacetime metric g^, background fields such as a Maxwellian-type 
gauge potential or a Kalb-Ramond gauge field (but not their gradients) 
and any relevant internal fields (that will be discussed in detail below). Note that 
the Goto-Nambu model has a constant Lagrangian density, namely 



-GN 



-Mo- (3.1) 



Upon infinitesimal variations in the background fields, and provided that indepen- 
dent internal fields are kept fixed (or alternatively that their dynamic equations of 
motion are satisfied), the action will change by 



5S 

where 



-\ J (T^Sg^ + W^SB^ + 2J»A li ) ^d 2 o , (3.2) 



T » v = 2 J^_ + Cr] »v (3.3) 



is the worldsheet the stress-energy tensor density, 



is the worldsheet electromagnetic current density, and 

w „u = 2 (3 5 ) 

OB^ 

is the worldsheet vorticity flux. The later will not be considered further in this 
paper. 

It useful to define two orthogonal unit vectors tangent to the worldsheet, one of 
them being time-like and the other space-like, such that 

- = V*% = 1 , = . (3.6) 

The eigenvalues of this ortho-normal frame are the energy density in the locally 
preferred string rest frame, which will henceforth be denoted by U, and the local 
string tension, denoted T (there should be no confusion with the vector defined 
in (3.6) and the stress-energy tensor T^ v defined in (3.3), respectively), so that one 
can write 

= UU^U" - TV^V" . (3.7) 
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Note that U and T are simply constants for a Goto-Nambu string, 



U = T = m , (3.8) 

but they are variable in general — hence the name 'elastic strings'. In particular, 
one should expect that the string tension in an elastic model will be reduced with 
respect to the Goto-Nambu case due to the mechanical effect of the current. 

Since elastic string models necessarily possess conserved currents, it is convenient 
to define a 'stream function' ijj on the worldsheet that will be constant along the 
current's flow lines. The part of the Lagrangian density C containing the internal 
fields is usually called the 'master function', and can be defined as a function of the 
magnitude of the gradient of this stream function, A = A(x), such that 

X = l ab Mb, (3.9) 

where the gauge covariant derivative is defined as 

= d a ^P - eA„x% . (3.10) 

Note that the definition of \ differs by a minus sign from that of Carter ||; the 
reason for this will become clear below. This 'dynamic' term contains charge cou- 
plings, whose relevance will be further discussed below. Nevertheless, whether or 
not these or other background gauge fields are present, it is always the form of the 
master function which determines the equation — or equations — of state. 

There is also a dual jlj| potential ■0, whose gradient is orthogonal to that of if), 
and the corresponding dual master function A = A(%) such that 

X = l a %ai>;b, (3.11) 

with the obvious definition for ip ;a . The duality between these descriptions means 
that the field equations for the stream function ip obtained with the master function 
A are the same as those for the dual potential ip obtained with the dual master 
function A. However, there will in general be two different equation of state relating 
the energy density U and the tension T; these correspond to what is known as the 
'magnetic' and 'electric' regimes, respectively corresponding to the cases 

Xmg < < Xmg (3.12) 

and 

Xel < < X el , (3.13) 

that are respectively characterised by space-like and time-like currents. In the 
degenerate null state limit, however, there will be a single equation, 

U = T = » . (3.14) 

Note that the distinction between a given model and its dual disappears in the 
absence of charge couplings; such models are then called 'self-dual' for obvious 
reasons. 

In each case the equation of state provides the expressions 

4 = | , (3.15) 

o dT v da , 

= (3 - 16) 
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for the extrinsic (that is transverse, or 'wiggle') and for the sound-type (longitudinal 
or 'woggle') perturbations of the worldsheet. Both of these must obey c 2 > (a 
requirement for local stability) and c 2 < 1 (a requirement for local causality) . These 
two speeds can be used to characterise the elastic model in question; in particular 
there is a straightforward but quite meaningful division of the models into supersonic 
(that is, those obeying ce > cl), transonic (eg = cl; only in the null limit is this 
common speed unity) and subsonic (eg < cl). 



B. Supersonic (superconducting) models 



Carter and Peter [|l0[ have recently proposed two supersonic elastic models to 
describe the behaviour of current-carrying cosmic strings. The Lagrangian density 
in the magnetic regime is 



A, 



-in 



2 , X 



X 



2fc n m 2 



(3.17) 



m a being the current carrier mass (which is at most of the order of the relevant 
Higgs mass); this is valid in the range 



1 

- - < 



X 



3 " k ml 
and obeys the equation of state 



< 1 



2m 2 + koml ' 



(3.18) 



U 



= 1 



kom 2 n 



fko~m a (T 



- 1 



fc m 2 



2 \ 1/2 



4m 2 y/2 m \m 2 8m 2 
On the other hand, the electric regime is described by the Lagrangian density 

fc m 2 



(3.19) 



A e ; = —m 



In 1 - 



which is valid in the range 



1 < 



X 



koml 



X 

koml J ' 



and the corresponding equation of state is 

U T koml \ ( 1 - Tim 2 
m z i \ koml/m^ 



(3.20) 



(3.21) 



(3.22) 



These models are supersonic for all space-like, and weak time-like currents, with 
the exception that in the null limit x — one has Cf, = ce = !• 



C. Equations of motion 



We now derive the microscopic equations of motion for elastic cosmic string mod- 
els. It is convenient to start by defining the quantity 

aa _ . 

Qab= A ab_ 2 ,;a,;b ,333) 
OX 
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then recalling the definition of x, (3.11), one can find the free string equations of 
motion in the usual (variational) way, obtaining 



(3.24) 



Also in a similar way to what was done in section [n], the effect of the frictional 
forces is accounted for by introducing a term ^—jF a on the right-hand side of 
(3.24). For exactly the same reasons as those of section, we will have 



pa 



(3.25) 



However, it should be remarked that the generalised definition of the friction length- 
scale for elastic models is 



i{ = 



A 



(3.26) 



(note that A is negative). 

Of course we now have a further equation for the scalar field tp, namely 



(3.27) 



Furthermore, the spacetime energy-momentum tensor and electromagnetic cur- 
rent will be given by 



and 



-gJ* 



-2c: 



- 7 9 a V>V (x - x(a, t)) d 2 a 



-lJ ab —^.aX%5(x-x(a,T))d 2 CT. 
ox 



(3.28) 



(3.29) 



The total string energy and charge in a spacetime where the l ine element is (2.19) 
are then defined as (we are again using the gauge choice ( 2.18 )) 



E = J d^v^T" = a j (-A+ tda 



and 



Q = J d 3 xV- (3 Vo = -2ea J ^~$fieda . 
The corresponding worldsheet charge and current densities are defined via 



and have the following values 



P = J 



-2e—eip ;0 
OX 



(3.30) 



(3.31) 



(3.32) 



(3.33) 
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Again, for the reasons explained above, a particularly relevant situation will be 
that of a chiral current, that is one in which 



This is equivalent to 



and therefore it implies that 



lab ff = • (3.35) 



ip' 2 = e 2 rjf , (3.36) 



dp = dj_ 
8t da 



(3.37) 



and that the total (spacetime) charge and current are also equal. Note that in the 
chiral case one also has 

dA 

X = 0, 2 0~ =1 > ( 3 - 38 ) 

so this is not equivalent to the Goto-Nambu case despite the fact that the equation 
of state is 

U = T = no- (3.39) 

Last but not least, one can always define the fundamental worldsheet current 
density vector 

i = = 2 v-7^7 W;b- (3.40) 

Although we have included the charge coupling term is the master function and 
its dual, it should be said that charge coupling effects are subdominant, and thus 
for most purposes they can be neglected (if nothing else, at least to a first-order 
approximation). This has been confirmed by Peter Jl6| , and is a consequence of 
the smallness of the coupling constants — for example, the electromagnetic coupling 
constant is e 2 ~ 1/137. In most of what follows we will therefore neglect the charge 
coupling. 

If an electromagnetic coupling does exist, it will be simply given by 

i" = , (3.41) 

where 

i" = x^ a , (3.42) 



is the corresponding spacetime current. Note that ( p.27[ ) is then just a Noether 
identity, 

i»=0. (3.43) 
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D. The chiral limit 



We now consider the (common) chiral limit of the two supersonic elastic models 
of Carter and Peter Effl, defined by the Lagrangian densities (3.17) and ( I3J20I ), 
respectively for the magnetic and electric regimes. Also, as we did in section O for 
the Witten model, we will interpret the charge coupling and the scalar field as being 
renormalised and neglect the coupling to external electromagnetic fields. 

Then, with our usual gauge choices and definitions of the d amping and friction 
length-scales, the microscopic string equations of motion (3.24) simplify to 



[e (1 + *)]•+ ^x 2 =*'-2-e*. 

t d a 



and 



e(l + *)x+— (l-x 2 )x = 
where ^ is defined as 



(I-*)- 



(3.44) 



* + 2-* x' + 2#x', (3.45) 
a 



Mo7oo 



(3.46) 



That is, these are exactly the same equations of motion as those of Witten's 
model (2.20-2.21) if one identifies the corresponding scalar fields, 



<f> = i/). 

Then, the worldsheet charge and current densities also coincide, 



(3.47) 



qeip, 



(3.48) 



(3.49) 



Finally, the total energy of a piece of string is given by 



E = fi a a / (1 + eda . 



(3.50) 



which we can immediately interpret as being split in an obvious way into a 'string' 
component and a 'current' component. This interpretation will be relevant below. 

Note that if we had preserved Carter's original sign conventions we would have 
found a difference of a factor of i between the two fields. But the important point is 
that the equality between the two theories in the chiral limit is not entirely trivial 
since, as we already pointed out, the motivations behind the build up of each of 
them are quite different. We have thus provided the first substantive evidence of 
the validity of the supersonic elastic models of Carter and Peter pi] . 



IV. CHIRAL LOOPS IN FLAT SPACETIME 



We will now study the evolution of current-carrying cosmic string loops, starting 
by considering the simplest case of circular loops in flat spacetime. We therefore 
choose the ansatz 
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x(t) = r (t) (sin 9, cos 9,0); 



(4.1) 



we also need an ansatz for the scalar field t/j (or <fi), which we will take to be 

i> = ^{T{t) + na)t c , (4.2) 

where the winding number per unit a, n, is a constant (due to the symmetry of 
our loop solution) and t c is a characteristic timescale — say the epoch of network 
formation. The chirality condition implies that 

eT = n. (4.3) 

Then the string equations of motion reduce to 

n 2 t 2 

= 1, (4.4) 



r+{l--fj r = 0, (4.5) 
together with the constraint 

M < \ ■ (4-6) 

Note that opposite signs of n correspond to left and right moving currents; naturally 
it always appears as n 2 in any relevant equation, and we will therefore be taking n 
to be positive. 

In figure [l] we plotted some relevant evolutionary properties of chiral supercon- 
ducting loops with different n's in flat spacetime. Note that these loops never 
collapse to zero size, and that their microscopic velocity is always less than unity 
(unlike in the Goto-Nambu case). Furthermore, there is a static solution with 

n = \ = T c f = 0; (47) 

in this case the energy is equally divided between the string and the current. 

It should also be noted that energy is transferred back and forth between the 
string and the current as the loop oscillates. We can easily determine the following 
quantities (the averages are over one oscillation period) 




(4.8) 



(4.9) 



1 



r- 



^ = -(l-4n 2 ), (4.10) 



2 

while the energy the string obeys 



;^™*> = i-„, (4.ii) 

-t^total 



12 



= 1-|»; (4.12) 

^total Z 

note that the energy of these loops is E to tai /t c = 2tt/j,. 

Finally, two other points that will have further relevance below. Firstly, a loops 
with a given conserved number n will reach a maximum microscopic velocity (and 
corresponding Lorentz factor) given by 

r 2 m ax = 1 - 4 ™ 2 . Imax = ^- (4.13) 

Secondly, for fixed n and initial velocity, there will be two possible choices of that 
can be made — the difference is that in one of them most of the energy will be in 
the string, while in the other it will be in the current. We will call these two cases 
the 'string branch' and the 'current branch'. In flat spacctime, the two choices give 
physically the same solution (they simply correspond to different initial phases of 
the oscillation), but this will not be true in general. 



V. CHIRAL LOOPS IN THE EXPANDING UNIVERSE 



The case of circular loops in expanding universes is analogous, and we keep the 
ansatze for x and tp, 

4> = ^{T{t) + na)t c . (5.1) 

The winding number per unit a and the function F are also constrained as before. 
In terms of these quantities the total energy of the loop can be written as 

Etotai = ^Q^total — M0 fl / ( 1 + j £C ^ 7 ^ ^-string + E current , (5-2) 

and the loops evolve according to 



r + (1 -r 2 ) 



r L ir I 2 . \ 1 



Id \ a 2 r 2 I r 



0. (5.3) 



It is convenient to define a macroscopic dimensionless parameter which, as we 
will show later, turns out to measure the loop's stability against collapse. We will 
define it by 

- 4wntc 

n=- ; (5.4) 

^total 

note that unlike n which is a constant for each loop, n is a variable parameter 
obeying 

0<n<l; (5.5) 

also n — corresponds to the Goto-Nambu case, while the n = 1 limit is the 
analogous of the flat spacetime static solution, here characterised by 

E s tring — E curren t , f = ; (^-^) 

in the approach to this limit one can easily establish that the loop's velocity (in the 
radiation epoch) and length in string evolve according to 
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V— = — ^string _ n . /r j\ 

t c 2 ' 2-Kt c x/] - 72 ' y ' ' 

these will be numerically confirmed below. 

An important difference with respect to the flat spacetime case is that now the 
string branch and the current branch (see figure || for a relevant particular case) 
represent two physically different solutions — something to be expected since damp- 
ing forces (that is, friction and expansion) act differently on the string and cur- 
rent energies. Since we will be mostly interested in chiral superconducting string 
loops formed in the friction-dominated regime (as no vortons will form in the 'free' 
regime) , we can safely assume that these loops are formed with zero velocity. Now, 
there is a very simple relation between Itotai , (string and n, namely 



(-total _ 2 
^-string Tl 



(l ± Vl-n 2 ^ . (5.8) 



The negative sign corresponds to the string branch, where as n goes from zero 
to unity we go from the Goto-Nambu case to the static case where the energy is 
split equally between the string and the current; the positive sign corresponds to 
the current branch, where the ratio of the energies in the strin g an d in the current 



decreases until it vanishes when n reaches zero again. Note that ( p.8|) can be inverted 
to give 



~X-strinq ( *~total -. \ ,~ 

n = 2- 9 - 1 . (5.9) 



1/2 

- 1 ' 

(total \ist ring 

In practice, it is not easily conceivable that in cosmological contexts loops can 
be formed with more energy in the current than in the string itself. Therefore, 
although for the sake of completeness we will be discussing the current branch in 
the remainder of this section, we will neglect it afterwards. 



Thus from (5.3) one obtains the evolution equation for £(£i,ti,rli,t), v((i,ti,rii,t) 
and other relevant quantities. As we will see below, a crucial quantity will be 
the the maximum velocity reached by each loop configuration during it evolution, 
v ma x((i, U,ni). If the loop does become a vorton, than its length will asymptotically 
be given by £ v (£i,rii). 

In figures ||-|^ we plot the cosmological evolution of some relevant GUT-scale 
chiral circular loops. We should mention that in order to save space, only one out 
of every forty points resulting from the numerical integrations is plotted, and this 
is the reason why some plots show irregularities. 

Figure || shows some relevant properties of the evolution of chiral circular GUT- 
scale loops formed at t = t c ; all have an initial total energy E to tai /27r/ii c = 10, but 
the distribution of the energy between the string and the current varies. 

Obviously, loops with higher currents will have smaller physical radii, and hence 
they will be less stretched by expansion and enter the horizon earlier, at which point 
they start oscillating — as can be confirmed in ||(a-b). Regarding the velocities, note 
the significant differences between loops in the 'string branch' (which still reach 
fairly high microscopic velocities, but never v — 1) and in the 'current branch' 
(which quickly become non-relativistic). Therefore the latter ones should definitely 
become vortons, and so it is perhaps fortunate that, as we pointed out above, we 
do not expect loops with such high currents to be produced in the early universe 
(at least, for GUT-scale networks). Note that in one of the cases shown the initial 
current is so high that the loop 'overshoots' and acquires a fairly large velocity, but 
friction quickly slows it down again. 

On the other hand, in the string branch the velocity is reduced with respect 
to the Goto-Nambu case, and a more detailed investigation will be needed to set 
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up some criterion defining which velocities will allow vorton formation — recall that 
relativistic velocities will imply charge losses and it will therefore be unrealistic to 
make any definite claims or predictions about such cases. 

The evolution of the fraction of the loop's energy in the current is particularly il- 
luminating (see ||(c)). This will obviously decrease while the loop is being stretched, 
and it will start oscillating when the loop falls in side the horizon. The oscillations 
are around the state with equipartition of the energy between the string and the 
current, which as we saw corresponds to a static solution in flat spacetime. Note 
that the effect of the friction force is to reduce the amplitude of these oscillations, 
so one can see that friction is in fact crucial for vorton formation. Naturally, loops 
with smaller velocities will undergo oscillations with smaller amplitudes, so again 
we confirm that these are the strongest vorton candidates. Finally, we have plotted 
the parameter n (which was defined in 5.4) in ||(d), and as one can easily see by 



comparison with the other three plots this is indeed a good indicator of whether 
or not a given loop can become a vorton — in fact, the 'phenomcnological' criterion 
that we mentioned above will be basically expressed in terms of the value of n onc e 
the loop is 'free' — that is, much smaller than the damping length defined in ( |2.22| ). 

On the other hand, radiative backreaction also tends to damp these energy os- 
cillations, and consequently increase n. Note that this has been shown to have the 
approximate form E = T em J 2 , and since T e m ~ 100 the timescale for this process 
is expected to be relatively short. 

Note that when loops become smaller than the damping lengthscale l& and reach 
the 'free' regime the following averages over a period hold (note that n becomes a 
constant in this limit — hence its usefulness) 

(f 2 ) = i(l-n 2 ), (5.10) 
<§^> = l-k (5.11) 

Etotal ^ 



F 2 1 

/ string \ <J 

^total * 



n; (5.12) 



the variance of the fraction of the energy in string is therefore 

A^^ = -n(l-n). (5.13) 

Etotal 4 

In figure [| we show chiral loops with the same initial conditions as |^, but starting 
to evolve at the epoch t+ when when friction becomes negligible . The differences 
are self-evident. Now, after a first period of growth of the total radius due to 
expansion, there is no mechanism forcing the loops to return this extra energy 
back to the medium when they fall inside the horizon. Consequently there is also 
no velocity damping (all loops will have microscopic velocities larger than 0.5) 
and the energy oscillations between the string and the current always have a large 
amplitude — so that n will never stabilise close to unity when the loops fall inside 
the horizon. 

Note that a loop with a very high initial current will again 'overshoots', but unlike 
in the case with friction here it can actually end up oscillating faster than another 
one in the 'current branch' but with a smaller current. This is because now there 
is no friction force that can damp this velocity overshoot. 

It can therefore be seen that vortons can only form during the friction-dominated 
epoch (as we expected) , and also that the earlier a loop is formed the larger will be 
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the region of the space of initial conditions that will originate them — because as we 
said the effect of friction is to increase n. Therefore, for cosmic strings formed at 
the GUT phase transition, the most favourable case for vorton formation is having 
the strings becoming superconducting at the GUT scale as well. We will use this 
assumption in the remainder of the paper. 

Finally, in figure || we plot the more realistic case of the evolution of GUT-scale 
loops having an initial string radius ten times smaller than the horizon, and different 
initial ratios of energies in the current and in the string — ranging from 10 -3 to 2. 

Now the total radius only suffers a small decrease, except in the case where 
one starts with n ~ 1, in which case velocity is so small that friction does not 
significantly affect the loop. Note that as n approaches unity we have v oc t~ l as we 
predicted, although for loops in the string branch there is an initial transient where 
v oc i -4 . Nevertheless, in the string branch loops do reach fairly high velocities 
during their first few oscillations, so that once more the issue of whether or not 
these become vortons is not entirely straightforward. 

Also note that for loops of this size the amplitude of the energy oscillations 
between the string and the current is negligibly small, except for the short transient 
period (typically lasting less than one Hubble time) for loops in the 'string branch' 
with fairly small currents. Clearly the relation between the initial conditions and 
the values of n and v needs to be looked at in more detail, and we shall do that in 
the next section. 



VI. CRITERIA FOR VORTON FORMATION 

In the previous section we saw that the evolution of chiral superconducting cosmic 
string loops depends sensitively on the conditions at formation. In particular, one 
would need to know in which cases one ends up with a vorton. 

Clearly, since we are not including radiative mechanisms at this stage, our crite- 
rion should be that loops whose velocity is always small (in a sense that will need to 
be made more precise) will become vortons, while those who are relativistic at some 
stage will suffer significant charge losses, so that their fate cannot be clearly asserted 
until a rigourous quantum-mechanical treatment of these processes is available. 

Thus we will explore in more detail the phase space of possible initial condi- 
tions in order to determine relevant properties of these loops. Figure |^ shows 
the maximum microscopic velocity v max (£i,ti,rii) reached by GUT loops formed 
at t = t c , 10 t c , 100 t c and i* ~ 855 t c , respectively; it is assumed that all such 
loops start their evolution with a negligibly small velocity — a reasonable assump- 
tion, since the network dynamics is friction-dominated until i*. In each case the 
horizontal axes correspond to the initial value of n and to the base-ten logarithm 
of the initial string radius relative to the horizon; recall that we only consider loops 
having initially most of their energy in the string (in other words, loops in the 
string branch). Note that the friction length-scale corresponds to about —1.5 in the 
vertical axis on the first plot, and to on the last (where it is equal to the horizon, 
by definition). 

It can be seen that any loop initially larger than the horizon will inevitably 
become relativistic. This is essentially because expansion will (temporarily, at least) 
decrease the fraction of the loop's energy in the current (and hence n). On the other 
hand, loops smaller than the friction length (and the horizon) have essentially no 
mechanism that can change n (neglecting radiation), so we will need fairly high 
initial currents in order to get non-relativistic velocities. 

Finally, for the case of loops being produced with sizes between the friction length- 
scale and the horizon, which is of course the cosmologically relevant case during the 
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friction-dominated epoch M , friction will force the loop to shrink (thereby increasing 
n), while the effect of the cosmological expansion will be small, so in order to have 
non-relativistic velocities we are allowed to have smaller initial values of n than in 
the previous case. 

From the analysis of figures ||-^| one can see that we need fairly high values of 
n when the loops reach the 'free' regime in order to have reasonable chances of 
producing GUT vortons in the 'string branch'. Now, according to Ref. the 
energy of a superconducting loop configuration with radius R is approximately 

W 2 

E = 2nfiR + 2-7rE— — = E str i ng + E current , (6.1) 
H 

1 12 

where W is the winding number and 2ttT,W ~ N L is the net particle number. The 
parameter E is the result of an integral over the string cross-section, it is a variable 
in general, but a constant in the chiral case, and expected to be of the order of the 
inverse of a coupling constant, £ > 20; we will in fact take E = 20 unless otherwise 
stated. This is minimised for a radius 



W \n 



(6.2) 



this corresponds to a vorton state. As expected, this minimum value is E v = 
2E string and corresponds to n = 1 . 

Now, suppose that the energy of a given configuration is a little higher than this 
minimum. That is, let E ~ (1 + x)E v . Then such a configuration will have a radius 

= 1 + ^ ± [(1 + x) 2 - 1] 1/2 ; (6.3) 

we will c hoo se the plus sign since it corresponds to the 'string branch'. Then we 



can use (5.9) to find the corresponding value of n: 



1 



. (6.4) 



1 + X + y/x(2 + x) 



These are useful expressions to introduce 'phenomenological' criteria for deciding 
which loop configurations will produce vortons. We note that these should be 
established on the basis of more detailed numerical studies of the microphysics of 
the currents; in particular, significant model-dependence is of course expected. 

As an example, if we take as a necessary condition for vorton formation that the 
energy of a given configuration is at most 10% higher than E v , we find that the 
value of n once the loop size becomes smaller than the damping length should obey 

Ufree > ~ 0.91 , (6.5) 

or equivalently that the average fraction of the loop's energy in the current must be 

/ E curren t , 



E, 



total 



) T >0.45; (6.6) 



Another (approximately) equivalent way of stating this is that a loop will not form 
a vorton state if it exceeds some maximum velocity v vor above which charge and 
current losses become effective. Note that a fast-moving loop will tend to develop 
cusps at which such losses should be particularly significant. Hence we require that 
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Vmax < Vvor if a given loop is to form a vorton. Of course v vor depends on x; for 
x = 0.1, we have 

v vor ~0.29. (6.7) 

Such a velocity limit is physically plausible, but a rigorous quantum mechanical 
treatment will be required to obtain more precise values. Note that the size of this 
vorton-forming region of parameter space is maximal at t c and decreases with time, 
vanishing not later than f*. 

If we choose less stringent criteria, say x = 0.5 or even x — 1.0, our bounds will 
respectively be 

2 

nfree > , , Vvor ~ 0.53 , (6.8) 



nfree > ^ , V vor ~ 0.61 ; (6.9) 



we will comment on the importance of the precise choice of x in section IX . 

Clearly, this only solves half of the problem — the other half is determining what 
exactly are the initial conditions at the formation of these loops, and in particular 
what are their currents. In other words, we need to know whereabouts in figure [j] 
do the loops form. This is a non-trivial problem, but we will discuss a simplified 
'toy model' for current evolution in the following section. 



VII. EVOLUTION OF THE CURRENTS 



Due to the strings' statistical nature, analytic evolution methods must be 'thermo- 
dynamic', that is one must describe the network by a small number of macroscopic 
(or 'averaged') quantities whose evolution equations are derived from the micro- 
scopic string equations of motion. The first such model providing a quantitative 
picture of the complete evolution of a string network (and the corresponding loop 
population) has been recently developed by the present authors JtJ, and we briefly 
summarise it here. 

We start by defining our averaged quantities, the energy of a piece of string, 

E = fia{r) J eda, (7.1) 

(e being the coordinate energy per unit a) , and the string RMS velocity, defined by 

o fx 2 edcr .„ . 

V =-J7^- (7 - 2) 

Distinguishing between long (or 'infinite') strings and loops, and knowing that 
the former should be Brownian we can define the long-string correlation length as 
Poo = fi/L 2 (see for an extensive discussion of these quantities, and others to 
be introduced below). A phenomenological term must then be included for the 
interchange of energy between long strings and loops. A 'loop chopping efficiency' 
parameter, expected to be slightly smaller than unity, is introduced to characterise 
loop production 

d Poo \ - Poo ,„ 

-7T ) =cv oc —. (7.3) 

at J to loops L 
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One can then derive the evolution equation for the correlation length L jTJ , which 
has the form 

2^ = 2HL(l + v 2 0O )+v 2 oo f-+cv oo ; (7.4) 

we point out again that the 'friction lengthscale' £{ will in general be that due to 
Everett scattering. 

One can also derive an evolution equation for the long string velocity with only 
a little more than Newton's second law 



— = 1 — v ) 
dt v ' 



*L- v (2H + ±. 



(7.5) 



here k is another phenomenological parameter that is equal to unity during the 
friction-dominated epoch and of order unity later . 

Finally, a careful analysis of the loop production mechanism leads to an expression 
for the energy density in loops. The idea is that at a given time one looks back at all 
the loops that have formed (and still have not decayed), finds their present lengths 
and then adds them together. Distinguishing between 'dynamical' and 'primordial' 
(that is, Vachaspati-Vilenkin) loops, we have 

Po{t) = J n dyn (t, t')£{t, t')dt' + j n pri (£', t)£ pri (£', t)dt . (7.6) 

Above £(t, t') is the length at time t of a loop produced at time t' (this will vanish 
if the loop has decayed), while ridyn(t,t')a 3 (t) = ni oop (t')a 3 (t'), where 

nioo P {t) = g^c^j (7.7) 

is the number of loops produced per unit time per unit volume. The factor g *~ 1 /V2 
accounts for the fact that not all of the energy lost by the long-string network ends 
up in the loops — part of it is lost by velocity redshift. We are assuming that loops 
produced at time t have an initial length £(t) = a(t)L(t) — in other words, that loop 
production is 'monochromatic' (see |7J for a discussion of this point). Similarly, for 
the Vachaspati-Vilenkin loops £ pr i(£' , t) is the length at time t of a loop formed with 
length £' , while n pr i(£' ,t)a 3 (t) = ni oop (£')a 3 (t c ), where ni oop (£') is the well-known 
Vachaspati-Vilenkin loop distribution. 

The above quantities are sufficient to quantitatively describe the large-scale char- 
acteristics of a cosmic string network. We will describe the evolution of the currents 
by a recently introduced toy model plf , which we now discuss in more detail. 

Our analysis will be based on the assumption that there is a 'superconducting 
correlation length', denoted £, which measures the scale over which one has coherent 
current and charge densities on the strings. Associated with this we can define N 
to be the number of uncorrelated current regions (in the long-string network) in a 
co-moving volume V = a 3 Lg as follows 



where L^, is the total long string length in the co-moving volume. 

Now, £ and N will obviously change in the course of the evolution of the string net- 
work, and we can immediately identify four possible sources of change — expansion, 
inter-commuting, loop production and internal dynamics on the string worldsheet. 
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We now consider each one of them. Firstly, we expect that in a co-moving volume 
the number of uncorrelated regions will not be affected by expansion, so 



dN 
~dt 



= 0. (7.9) 

expansion 







Now consider the effect of intcr-commutings (whether or not a loop is produced) . La- 
guna and Matzner [ [17] have numerically shown that whenever two current-carrying 
strings cross, they inter-commute and a region of intermediate current is created. 
This means that inter-commutings will in general create four new regions (see figure 
(a)). Since, according to our analytic evolution model the inter-commuting rate 
is 

^L) = l?Z°L (7 10) 

/ inter commuting 

we immediately obtain the following effect on N 



/ inter commuting 



again this assumes that loops have a size £(t) = a(t)L(t) at formation, and that 
once the long-string network reaches the linear scaling regime we have a(t) — a sc = 
const, (sec j7J). 

However, an important correction is necessary to account for the fact that when 
regions with size of order £ or smaller self- intersect it is possi ble (s ee figure (b-c)) 



that no new regions are produced. Thus we must multiply (7.11) by a correction 
factor 



Fxll) =1,§>1 (7.12) 



■2c [ i - i 



i 



(l_ a) _|<l. (7.13) 



The slightly complicated behaviour of F\ is nevertheless easy to understand. 
The point is that numerical simulations show that there are two types of inter- 
commutings. Firstly, 'large-scale' ones always occur at a scale L; a fraction a of 
the inter-commutings should be of this type. If this happens between two long- 
strings (that is, no loop is produced) we always expect to create new regions, since 
there is no reason for currents in different 'infinite' stings to be correlated. On the 
other hand, if what we have is a long string self-intersecting to produce a loop of 
size smaller than £ (a fraction 2c of these inter-commutings should produce loops), 
then we might not form new regions — for each length, the fraction of these self- 
intersections that produce new regions is essentially given by the ratio of the size 
of the region and the superconducting correlation length. The remainder of the 
inter-commutings are associated with the presence of small-scale structure on the 
strings, and occur by repeated self-intersections of a given string so the £/£ cut- 
off always applies. Notice that the second term vanishes if a = 1 (as it should) 
but it rapidly becomes dominant as a starts deviating from unity. Also note that 
the overall inter-commuting effect is approximately a-independent (more on this 
below) . 

Of course, when the inter-commuting does produce a loop, the regions in the 
corresponding segment are removed from the network, together with one of the 
newly created 'intermediate' regions, and we similarly have 
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where the analogous correction factor F2 is of the form 



' (7) 1 ( 7 - 15 ) 



Fi[ 7 ) = 7 ,f <1. (7.16) 



Note that string length is always removed from the long string network when loops 
form, regardless of whether or not current regions are. This is in fact the main effect 
of loop production, as can be seen by noting that ( [7.14 ) is approximately indepen- 



dent of the parameter characterising the loop size, a. In the friction-dominated 
regime, a is of order unity, and when it becomes much smaller (in the free regime) 
the a-dependencies in the numerator and in the denominator cancel out. One can 
readily see that this is physically plausible: when a ~ 1 (in the friction-dominated 
epoch) few loops are produced, but each one of them removes a significant number 
of regions; on the other hand, when a is small, many more loops are produced, but 
only a few of them will remove regions. 

Finally, there is the dynamic term. When regions with opposite currents inter- 
commute, new charged regions are created, setting up alternate currents. One ex- 
pects electromagnetic processes to make these currents die down, so that the charged 
region will eventually equilibrate with its neighbours. The simulations of Laguna 
and Matzner Jl7|] provide qualitative support for this intuitive picture. Clearly, 
this indicates that some kind of 'equilibration' process is effectively acting between 
neighbouring current regions, which will counteract the creation of new regions by 
inter-commuting. While it is beyond our means to derive an 'equilibration term' 
from first principles we will, as a first approximation, introduce a phenomenological 
term. We will model this current decay by assuming that after each Hubble time, 
a fraction / of the N regions existing at its start will have equilibrated with one of 
its neighbours, 



~dt) 



) = —fHN; (7.17) 

dynamics 



note that new regions are obviously created by inter-commuting during the Hubble 
time in question, so that / can be larger than unity. Alternatively we can say 
that for a given /, the number of regions in a given volume at a time t will have 
disappeared due to equilibration at a time t + (/if) -1 . We therefore obtain the 
following evolution equation for N 



where we have re-defined the correction factor 



a[i)=2-c(i + 2),l > 1 (7.1!)) 



G f -J = 2(1 - 2c)a + (2 - 3c - 2a + 4ac)-, | < 1 . (7.20) 

Note that when £ 3> ^ the net effect of inter-commuting and loop production is to 
remove uncorrelated regions (because each loop formed removes a large number of 
them); otherwise, the net effect is to create new regions. 
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However, for what follows it is convenient to re-write it in two alternative forms. 
Firstly, we can define Nl to be the number of uncorrelated current regions per 
long-string correlation length, 

Nl = f = N TZ ; (7 - 21) 

this is useful because, as was first pointed out by Davis and Shellard ||, we expect 
the net charge of a superconducting loop to be given by 

Q ~ eN x L /2 . (7.22) 



In terms of Nl, (7.18) has the form 

X L (3t& f)HN L + l V ^N L +( -G(aN L ) + ^N L ) : ! 7.2:; ) 



dt v 00 " 2 ii \a y ' 2 

note that to obtain this one ne eds to substitute the evolution equation for the 
long-string correlation length L (TA), and that one can equivalently define G as 

G{aN L ) = 2-S(aN L + 2),aN L > 1 (7.24) 
G (aN L ) = 2(1 - 2c)a + (2 - 3c - 2a + 4ad)aN L ,aN L < 1 . (7.25) 

Yet another useful form follows from defining Nh to be the number of uncorre- 
lated current regions in one Hubble volume, 

*H = If = NL%i (7-26) 



in this case we have 



^ = (3 - f)HN H + G(aN L ) V ^. (7.27) 
dt a Li 1 



VIII. THE IMPORTANCE OF EQUILIBRATION 

Now the question is, of course, what is /? From a more intuitive point of view, 
an equivalent question is the following: given a particular piece of string with a 
given current, is it more likely to disappear from the network by this equilibration 
mechanism or by being incorporated in a loop? Even though a precise answer 
can probably only be given by means of a numerical simulation, some very simple 
physical arguments can be used to restrict it. We should point out, however, that 
many of the results of the following sections do not depend crucially on the value 
of/.' 

Firstly, correlations cannot obviously be established faster than the speed of light 
(that is, we must have £ < t), so that we should impose that 

this leads to an upper bound on /, which we can write, defining L = jt, as 

fmax(c, 7, Voo) =3 + 4(1-25)^. (8.2) 

7 
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(In this section we will concentrate on the bounds on / in the radiation epoch — 
analogous results can obtained for the matter epoch.) We explicitly write the de- 
pendencies of fmax to emphasise that this is the maximum value of / which satisfies 



1) for a given set of properties of the cosmic string network. 



On the other hand, if an equilibration mechanism such as that modelled by (7.17) 
exists [ H , it is reasonable to assume that it will prevent Njj from growing without 
limit — possibly through a backreaction mechanism as in the case of gravitational 
radiation for wiggly Goto-Nambu strings — and eventually it will make it become 
constant (meaning that £ is scaling linearly). In other words, we can assume that 
there should be a large Njj (which we need not specify) such that 



<0; (8.3) 



we can therefore find a lower bound on / which satisfies this, 

v 

fmin{c, 7, Voo) = 3 - 2c— . (8.4) 

7 

Again this varies as the n etwork evolves. Note that the crucial point about this 



construction is that ( 7.17 ) depends linearly on Nh- 

Using the quantitative evolution model of the present authors we have plotted 
fmax and fmin during the friction-dominated epoch with initial conditions typical 
of first and second order phase transitions, in figure ^. These plots are fairly easy 
to interpret. Perhaps the most surprising result is the large values of f m ax allowed 
when the long-string correlation length is well below the horizon. This is because 
in this case the loops chopped off by the network are small, so that each one of 
them removes relatively few current regions — one could therefore have an extremely 



efficient equilibration mechanism and still obey the constraint (3.1). Hence we can 
see from figure |^ that if our toy model, and in particular the ansatz (7.17) is valid 
the constraints on / are much stronger for a first-order phase transition. One can 
also see that / = 3 is the only value that is acceptable at all times, regardless of 
the initial condition. However, it is at present unclear if there is something 'special' 
about this value. A numerical simulation is presumably the only way to clarify this 
issue. 

Note that any value f m %n < / < fmax once the network reaches the linear scaling 
regime 

1.88 ~ 3 - -^L < / < 3 + - fj. 1 , ~ 2 f ^ ~ 22.4 , (8.5) 
k + c k 1 / 2 (k + c) 3 ^ 2 

leads to a constant value of Nl and that this corresponds to £ scaling as the long- 
string correlation length L. Different values of / lead to different scaling values of 
Nl (with larger /'s corresponding to smaller Nl's as expected) at least in some 
region of the space of initial conditions, but for any / ^ 3 one can think of some set 
of physically viable initial conditions for which either causality would be violated 
at some stage of the evolution or the number of uncorrelated regions would grow 
without bound. The scaling value of Nl can be written in terms of the properties 
of the string network as 

one can see that in this regime the / dependence is rather weak, unless / is just 
above f min . 
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We emphasise that while the f m ax bound is unavoidable (being a consequence of 
causality), / m j„ is less robust and could well be disproved by a detailed numerical 
study. Therefore, in what follows we will discuss two cases, / = and / = 3 
which should represent the scenarios of ineffective (or non-existent) and effective 
equilibration. 

For a given /, we can now solve (7.23) numerically, coupled with the evolution 
equations for the long-string correlation length and average velocity (see |7j]). This 
therefore allow s us t o know the size of the loops formed by the network at each time 
and (through ( [7.22 )) the initial current they will carry. On is then in a position 
of applying the criteria established in section VI in order to decide whether or not 
each loop will form a vorton. 

We should also say at this stage that once the network leaves the friction- 
dominated regime and strings become relativistic other mechanisms (notably ra- 
diation) can cause charge losses in the long strings (as well as in loops). Hence our 
toy model can at best provide order-of-magnitude estimates in this regime. On the 
other hand, we expect it to be quite accurate (pending a more detailed numeri- 
cal study) in the friction-dominated epoch — which is of course relevant for vorton 
formation. 



IX. GUT-SCALE ANALYSIS 



In figures pf|lC| we plot the result of the numerical integration of (7.23), for initial 



conditions representative of string-forming and superconducting phase transitions 
of first and second order, for the cases / = and / = 3. We are assuming that these 
occur at around the same (GUT) energy scale since, as we have shown in section 



VI , this is most favourable situation for vorton formation. It was also assumed that 
the value of a in the linear scaling regime is a sc ~ 1CT 3 (see ref. [§). 

The differences between the two cases are considerable. Firstly, if there is no 
equilibration mechanism (/ = 0, see figure ^), the number of uncorrelated regions 
per long-string correlation length Nl, never decreases. In this case there are simple 
scaling laws for Nl and £. One finds that £ is conformally stretched during the 
stretching regime (just like the long-string correlation length, L oc t 1 / 2 ), and so Nl 
is approximately constant. However, as inter-commutings start creating new regions 
Nl begins to increase, and it grows as t 3 / 2 during the Kibble regime (where L oc i 5 / 4 , 
so £ tx t" 1 / 4 ). Finally, once the network reaches the linear scaling regime, L cx t, 
the number of uncorrelated current regions grows as Nl oc t, which corresponds to 
£ oc const. As expected, in this case the network keeps a 'memory' of its initial 
conditions. 

On the other hand, if there is an efficient enough equilibration mechanism (see 



figure f for the case / = 3) then Nl decreases while the network is being confor- 
mally stretched. In the Kibble regime, the increased number of inter-commutings 
again drives Nl up, and after a has evolved into its linear regime value, £ itself 
reaches a scaling value and hence Nl becomes a constant. In the intermediate case 
of a small but non-zero /, Nl decreases during the stretching regime but grows 
without limit afterwards, and the precise values of the scaling laws depend on /. 
Also, the network will preserve a 'memory' of the order of the string-forming phase 
transition, but not of the order of the superconducting one. 

This therefore solves the other half of our problem. Knowing the loop size at 
formation at all times R| at the typical current that each loop carries at that epoch 
(from the above toy model) one can then apply some criterion (possibly of the 



type discussed in section VI ) to decide which loops have a reasonable possibility of 
becoming vortons. 
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Our quantitative string evolution model Q| allows us to determine the size o f th e 
loops formed at each epoch, l(t) = a(t)L(t). On the other hand, according to ( |5.9|) , 
to find the initial n we need to know the ratio of the energies in the string and in 
the current. Now, the energy of a superconducting loop configuration with a radius 
R is given by fl6.1| ), so after some algebra we find 

E current _ 167T 3 7V" G[lN L 

E stHng ~ 45S a 2 j 2 x 2 ' 1 ' > 

where L = jt, t = xt c and AT is the number of effectively massless degrees of free- 
dom. Note that the minimum value of 7 is of the order of (G/Lt) 1 / 2 (but slightly 
larger — see (7)), so the crucial factor in this equation, and hence for vorton forma- 
tion, is how much Nl can grow. This alone tells us that the higher the energy scale 
at which the string network forms, the less likely it is to produce vortons, since it 
will be friction-dominated (and hence non-relativistic) for a shorter period of time. 
In order to make vortons, loops should be formed with a high enough Nl to allow 
them to remain non-relativistic thereafter — otherwise, they will eventually become 
relativistic and hence liable to charge losses. As we already pointed out, making the 
strings become superconducting sometime after they form does not help — it merely 
reduces the time available to build up charges and currents. 

Contrary to current belief (which is based on rather more qualitative estimates) 
we do not expect any vortons to be produced by GUT-scale cosmic string networks. 
In figures |ll|-|l2| we plot the paths of initial conditions in n-R space for dynamic 
and Vachaspati-Vilenkin loops Jig] formed during the friction-dominated epoch in 
the cases / = and / = 3. For the 'dynamic' loops, we only plot loops formed until 
100 t c (notice that n decreases after this epoch). We consider initial conditions for 
the string network that are characteristic of first- and second-order string-forming 
and superconducting phase transitions. Note that the difference in the initial n's 
between the two cases is smaller than the difference between the corresponding 
Nl's; this is because n is approximately proportional to N L . 

One can see that, even if we choose the less stringent of our three suggested 
criteria, calling a vorton any loop configuration with an energy up to twice the 
minimum value (that is, x = 1) we still get no GUT vortons. In fact, one would 
need to choose a limiting velocity v vor ~ 0.7 for GUT vorton production to occur 
in this model — and even so, only in the case when equilibration is efficient and the 
string-forming phase transition is of second order. However, we should emphasise 
the issues of the precise vorton formation criterion, as well as that of the value 
of the 'equilibration parameter' /, can only be settled by means of more detailed 
numerical studies of the microphysics of these loop configurations. 

This is an appropriate point at which to add a cautionary note about the 
quantum-mechanical stability of vortons. This is a rather involved and model- 
dependent question which has been briefly discussed in ref. 0. The vorton gains 
additional stability because its charge carriers must tunnel off the string by tak- 
ing both charge and angular momentum; the larger a vorton is, the more stable it 
is. With electromagnetic fields present, pair creation provides an alternative decay 
mechanism, but for a chiral vorton with null fields (E 2 — B 2 ) near the string, this 
mechanism is strongly suppressed. As the chiral state tends to be an attractor for a 
wide range of initial conditions || this again encourages us to believe that vortons 
should generically be quantum mechanically stable. Nevertheless, one can make 
special parameter choices for which vorton lifetimes are very brief, notably when 
the string and current-forming phase transitions are widely separated in energy 
scale. This subject clearly deserves a more thorough investigation. 
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X. CALCULATING VORTON DENSITIES 



Vorton densities can be calculated using a fairly straightforward modification to 
the method developed in M and summarised above. We now have 



The original model for Goto-Nambu strings included an averaged evolution equa- 
tion for the length I of each loop, which made the above calculation relatively 
easy. Here, an analogous averaged equation for a superconducting loop is presently 
unavailable, but the loop size (and velocity) can be determined by evolving the 
microscopic equation of motion (^^). The functions Wi(t') and W2(i') are 'win- 
dow' functions — typically combinations of Heaviside functions — selecting the time 
interval in the evolution of the network (and the interval in the length of Vachaspati- 
Vilenkin loops) which will produce vortons, according to the particular criterion that 
one chooses to impose. Notice that these will depend on a number of parameters, 
including the initial conditions of the cosmic string network. Also, they should in 
principle include a factor accounting for the fact that it takes some time for each 
loop to reach a vorton configuration (that is, even if a given loop will eventually 
form a vorton, it should not be included in the vorton density until some time after 
it is 'chopped off' from the long-string network). However, note that figure ^| seems 
to indicate that this evolution, if it happens at all, is quite fast — it takes less than 
a Hubble time. 

Note that although in the evolution of the loops the effects of the currents are 
properly accounted for (with the exception of radiative mechanisms), the evolution 
of the long string network doesn't take into account of possible effects of the build- 
up of the currents. Still, we expect the neglect of these effects to be a reasonable 
assumption. This is because such effects should only become important (if ever) at 
late times when the network has had time to build up large currents while, as we will 
shortly see, most of the energy density in vortons is produced fairly soon after the 
network forms (but a possible exception to this can occur if there are background 
magnetic fields which can increase the current build-up rate). 

Thus calculation of vorton densities is a two-stage process. Firstly, one must 
study the microphysics of the particular model that one is interested in, in order to 
derive its microscopic equations of motion and in particular to construct appropriate 
expressions for the 'window functions' W\{t') and W^') which will determine at 
which stages of the evolution of the string network one can form vortons. Secondly, 
one can use the velocity-dependent one-scale model and the model for the evolution 
of the currents on the long strings, together with the microscopic loop equations 
of motion (or an averaged version of them) to determine the vorton density using 
Ql0.2j ). Typically there will be a single time interval t c < t star t < t < t stop < 
at which vortons will form, but it is relatively easy to think of initial conditions 
for which vortons form at two different time intervals. Values of t star t and t stop in 
specific models will be discussed in a forthcoming publication ]l2|] . 

Since, as we pointed out, there is some uncertainty in some crucial parameters 
of this model, we will limit ourselves in this paper to calculate the vorton density 
for GUT and electroweak string networks in the 'best' (or 'worst' according to 
opinion) possible case where there is no equilibration (that is, / = 0), the string- 
forming and superconducting phase transitions are both of second order, and all the 
loops produce vortons (hence our criterion is simply n/ ree > 0). Notice that this 




(10.1) 



(10.2) 
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last condition is unrealistic for GUT networks (where, as we already indicated, we 
don't expect vortons to from) but is plausible for electroweak networks. Still, we 
will assume that vortons can only be formed while the network in in the friction- 
dominated regime. Also, since one presumably needs to have quite efficient radiation 
mechanisms for all loops to relax into vortons, we will assume that such relaxation 
is instantaneous — thus W% = 1, while W\ is unity in the friction-dominated epoch 
and vanishes afterwards. 

Figure |l3| displays the resulting vorton densities, relative to the background and 
matter densities. Firstly, we confirm that most of the energy density in vortons is 
produced soon after the network forms. In the case of GUT strings, we see that 
vortons would only dominate the energy density of the universe about four orders of 
magnitude in time after the epoch of network formation, that is soon after friction- 
domination ends (recall that for a GUT network t+ ~ 850i c ). Thus even if all 
these vortons formed, they wouldn't contradict the standard cosmological scenario 
provided that they decayed soon after t*, when the network becomes free. In any 
case we emphasise that this 'worst case' scenario is not realistic for GUT-scale 
strings, and indeed (as discussed previously) we do not expect GUT-scale vortons 
to form at all. 

On the other hand, electroweak string networks are friction-dominated until after 
the radiation-matter transition, so the vorton density has been slowly building up 
relative to that of matter until very recently. We find that this density today would 
be about 6% of the critical density. On the other hand, a string network formed 
at T ~ 10 4 GeV would provide a maximal vorton density equal to the critical 
density. This is therefore the strongest possible vorton constraint — it is based on 
the assumption that all loops form vortons. Naturally, realistic models are not 
expected to be fully efficient in producing vortons, and furthermore the relevant 
phase transitions are not necessarily of second order. One can therefore conjecture 
that the dark matter problem might be solved by a superconducting string network 
formed at an energy scale of T ~ 10 5 — 10 6 GeV. Note that there are a number of 
super-symmetric models producing such networks (see for example P90. We will 
present a more detailed analysis of these issues in a future publication JL2] . 



XI. CONCLUSIONS 



In this paper we have presented the first rigourous study of the cosmological 
evolution of superconducting strings in the limit of chiral currents. We have shown 
that in this limit the elastic string model of Carter & Peter Jl(| coincides with the 
model derived from first principles by Witten [jjj . 

By analysing physically relevant loop solutions of the microscopic equations of 
motion for these strings, we have verified that the effect of frictional damping is 
crucial for vorton formation. We then defined suitable parameters characterising 
the evolution of these loops, and in particular whether or not they become vortons. 
In particular, we have established the usefulness of the 'stability parameter' n. In 
general, it is more difficult to form vortons when the string-forming phase transition 
is of first order. This is because such networks produce, during their evolution in the 
stretching regime, loops with a size close to that of the horizon; these will therefore 
be significantly affected by expansion, which tends to decrease the fraction of the 
loops's energy in the current — whereas friction tends to increase it. 

After introducing a simple 'toy model' for the evolution of currents on the strings 
p"l| ] , we have considered the cases of first and second-order GUT-scale string- forming 
and superconducting phase transitions (which is the most favourable GUT case of 
vorton formation since frictional forces can act longer). We have presented evidence 
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suggesting that GUT-scale string networks might well produce no vortons, and that 
even if they do, this will not necessarily rule out such models. This is in contradic- 
tion with previous, less detailed studies and hence calls for a re-examination 
of a number of cosmological scenarios involving superconducting strings. Notably, 
these strings could be at the origin of the observed galactic magnetic fields . 

Finally, we have explicitly calculated the vorton density in two 'extreme' cases to 
illustrate the method that one should follow once the microphysical properties of 
these networks are known in more detail. For electroweak-scale string networks, we 
have found that vortons can produce up to about 6% of the critical density of the 
universe. On the other hand, it is conceivable that superconducting string networks 
formed at an energy scale T ~ 10 4 — 10 6 GeV (depending on details of the model) 
can solve the dark matter problem. 

The detailed analysis presented in this paper for GUT stings can obviously be 
extended to other energy scales — this will be the subject of a forthcoming publica- 
tion Jl^] . Obviously, as we lower the energy scale, the frictional force becomes more 
and more important and acts for a longer time. Hence the vorton-forming region of 
parameter space increases, and by the electroweak scale almost all loops chopped 
off the long-string network will become vortons. We therefore conclude that in ad- 
dition to the low-G/i regime (which as we saw includes the electroweak scale) where 
vortons can be a source of dark matter and to an intermediate-G^ range in which 
vortons would be too massive to be compatible with standard cosmology (thereby 
excluding these models) , there is also a high-G/z regime (of which the GUT scale is 
part) in which vortons don't form at all and therefore no cosmological constraints 
based on them can be set. It is then curious (to say the least) that vorton con- 
straints can be used to rule out cosmic string models in a wide range of energy 
scales G/i, but not those formed around the GUT or the electroweak scales, where 
cosmic strings can be cosmologically useful. 
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FIG. 1. The flat spacetime evo luti on of chiral circular string loops characterised by a 
conserved quantity n (defined in 4.2) having the value (solid lines), 0.1 (dashed), 0.3 
(dot-dashed), 0.45 (dotted) and 0.5 (star). Note that the first corresponds to a simple 
Goto-Nambu lo op, while the last is a static solution. Plots respectively show r as function 
of r (defined in 4.1) (a) and the fraction of the energy in the string as a function of r (b). 




FIG. 2. The logarithm of the ratio of the energies in the current and in the string for 
chiral circular loops with zero velocity, as a function of the logarithm of parameter n. Note 
that there are two different branches, hereafter called the 'current branch' (top) and the 
'string branch' (bottom). 
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FIG. 3. The evolution of chiral circular GUT-scale string loops formed at t — t c . 
All loops have an initial total energy Etot/2irfJ>t c = 10, but different initial string en- 
ergies — respectively 9.3, 7.2, 5.0, 2.8 and 0.7; in (a) and (d), these are respectively shown 
in solid, dashed, dash-dotted, dotted and starred lines, while in (b) and (c) they are cor- 
respondingly shown by lighter shades of gray. Plots show the total energy Etotai relative 
to t c (a), the microscopic velocity (b), the fraction of the loop's energy in the current (c) 
and the parameter n defined in 5.4 (d). 
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FIG. 4. The evolution of chiral circular GUT-scale string loops formed at t = t*. 
All loops have an initial total energy Etot/2nfJ,t* = 10, but different initial string en- 
ergies — respectively 9.3, 7.2, 5.0, 2.8 and 0.7; in (a) and (d), these are respectively shown 
in solid, dashed, dash-dotted, dotted and starred lines, while in (b) and (c) they are cor- 
respondingly shown by lighter shades of gray. Plots show the total energy e to tai relative 
to t* (a), the microscopic velocity (b), the fraction of the loop's energy in the current (c) 
and the parameter n defined in 5.4 (d). 
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FIG. 5. The evolution of chiral circular GUT-scale string loops formed at t = t c . All 
loops have an initial string energy E a trin g /2-Kfit c = 0.1, but different initial ratios of ener- 
gies in the current and the string — the cases 1CT 3 , 1CT 2 , HT 1 , 1.0 and 2.0 are respectively 
shown in solid, dashed, dash-dotted, dotted and starred lines. Plots show the total energy 
Etotai relative to t c (a), the (base-ten) logarithm of the microscopic velocity (b), the frac- 
tion of the loop's energy in the current (c) and the (base-ten) logarithm of the parameter 
l-n(d). 
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FIG. 6. The maximum microscopic velocity reached by circular GUT-scale chiral super- 
conducting string loops formed at a time tf = t c ~ 10~ 39 s (top left) and tf = 10, 100, 855 t c 
(clockwise), the later being t* The 'x' axis corresponds to the initial value of the parameter 
n, going from zero (the Goto-Nambu case) to unity; in the 'y' axis the base-ten log of the 
string radius relative to the horizon size goes from —2 to 1. Note that in the first graph 
the friction lengthscale corresponds to log£f/t ~ —1.5, while in the last one log If /t = 0. 
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FIG. 7. Some relevant inter-commuting configurations. The arrows mark the limits of 
regions with correlated currents. Plot (a) shows a typical inter-commuting creating four 
new current regions, while (b-c) show than on scales smaller than the current correlation 
length loop production may (c) or may not (b) remove current regions from the long-string 
network. 
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FIG. 8. The evolution of / m i„ (lower pair of curves) and fmax (upper pair of curves) for 
first order (solid lines) and second order (dotted lines) string-forming phase transitions. 
Time is in orders of magnitude from the epoch of string formation. 
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FIG. 9. The evolution of the number of uncorrelated current regions per long-string 
correlation length, Nl, for the case / = (the bottom plot is a friction-dominated epoch 
close-up of the top one) assuming that the orders of the string-forming and superconducting 
phase transitions are respectively: 1st & 1st (solid lines), 1st & 2nd (dashed), 2nd & 1st 
(dash-dotted) and 2nd & 2nd (dotted). Time is in orders of magnitude from the epoch of 
string formation. 
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FIG. 10. The evolution of the number of uncorrelated current regions per long-string 
correlation length, Nl, for the case / = 3 (the bottom plot is a friction-dominated epoch 
close-up of the top one) assuming that the orders of the string-forming and superconducting 
phase transitions are respectively: 1st & 1st (solid lines), 1st & 2nd (dashed), 2nd & 1st 
(dash-dotted) and 2nd & 2nd (dotted). Time is in orders of magnitude from the epoch of 
string formation. 
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Dynamic loops, f=0 
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FIG. 11. The initial conditions for loop formation in ri-R space, for the case / = 0, 
assuming that the orders of the string- forming and superconducting phase transitions are 
respectively: 1st & 1st (solid lines), 1st & 2nd (dashed), 2nd & 1st (dash-dotted) and 2nd 
& 2nd (dotted). The top plot corresponds to dynamic loops formed between t c and 100 t c 
(in the first two curves loops formed at t c are at large R; in the later two they are at small 
R). The bottom corresponds to 'primordial' Vachaspati-Vilenkin loops formed at t c and 
having lengths between L c and 10 L c . 
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Dynamic loops, f=3 
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FIG. 12. The initial conditions for loop formation in ri-R space, for the case / = 3, 
assuming that the orders of the string- forming and superconducting phase transitions are 
respectively: 1st & 1st (solid lines), 1st & 2nd (dashed), 2nd & 1st (dash-dotted) and 2nd 
& 2nd (dotted). The top plot corresponds to dynamic loops formed between t c and 100 t c 
(in the first two curves loops formed at t c are at large R; in the later two they are at small 
R). The bottom corresponds to 'primordial' Vachaspati-Vilenkin loops formed at t c and 
having lengths between L c and 10 L c . 
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FIG. 13. The maximum possible vorton densities relative to the background (solid lines) 
and ordinary matter (dotted lines) densities, for GUT and electroweak-scale string net- 
works. Time is in orders of magnitude from the epoch of string formation; the plots end 
at the present epoch. 
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